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1.  Introduction 

This  note  is  partially  motivated  by  a paper  by  Krakov/ski  (1973)  in  which 
he  describes  the  following  process.  At  time  zero  an  individual  is  diosen  at 
random  from  a -oopulation  A of  new  (age  zero)  individuals.  'n'.en  this  indivi- 
dual fails,  it  is  replaced  by  an  individual  chosen  at  random  from  the  survivors 
among  a population  B created  (i.e.  at  age  zero)  at  t’le  same  time  as  A. 

Kralcowski  found  an  expression  for  the  distribution  of  failure  time  of  this 
replacement. 

!'fe  have  developed  this  model  to  allo!^'  for  Kxjre  tJian  one  stage  of  replace- 

j 

ment,  v/ith  different  nopulations  providing  the  replacement  at  different  stages.  ' 

Tlie  condition  tiiat  replacements  are  diosen  from  populations  aging  simul- 
taneously is  retained.  ! 

Furthermore,  (moti'/ated  by  B.C.  Arnold's  (1975)  v/ork  on  the  classifica- 
tion of  multivariate  exponential  distributions  based  on  hierarchal  successive 
damages  v/hich  extends  the  v/ell  loio'.n  fatal  shock  model  for  iiarshall  and  j 

I 

Dikin' s (1967)  bivariate  exponential  distribution)  we  have  imposed  a hierarclial  j 

structure  on  the  replacement  scheme  in  tvra  different  ways,  corresponding  to  j 

tit)  different  replacement  strategies  for  a nur-ber  of  elements  which  are  in  | 

I 

service  simultaneously.  ] 


it 
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In  all  these  cases  we  are  primarily  interested  in  the  final  failure 
time,  but  we  also  derive  the  joint  distribution  of  failure  (and  replacement) 
times  at  intermediate  stages. 

2.  Hierarchal  Replacement  Systems  (Procedures) 

We  will  nov;  describe  the  two  replacement  systems  (procedures)  which  we 
will  study,  first  descriptively  and  then  more  formally  introducing  a special 
notation. 


System  (Procedure)  1:  We  suppose  that  there  are  v ‘'positions"  each  occupied 
by  an  "element"  from  a given  set  A.  Initially  all  k positions  are  serviced 
by  a single  "component."  The  exponents  are  individuals  which  are  assumed 
to  be  cliosen  at  random  from  a population  with  survival  distribution  function 
(SDF)  SQ(t).  When  a certain  component  fails,  it  is  replaced  by  Sj^  conponents 
each  servicing  a disjoint  subset  of  kj^,...,kg  of  the  k positions  (kj^+k2+... 
■•■ke  = k) . Tnese  new  components  are  chosen  at  random  from  a population  which 
vreis  originated  ("bom")  at  the  same  time  as  the  original  component,  but  has 
survival  distribution  Sj^(t). 

Each  of  the  Sj^  subsets  nov/  has  the  same  status  as  the  original  set  of  k 

2 

elements.  The  i-th  subset  is  divided  into  s.^  sub- sets,  and  v/hen  the  com- 

i2  ’ 

ponent  servicing  this  subset  fails  it  is  replaced  by  s.2  components  chosen 
at  random  from  a population  with  initial  survival  distribution  function  (SDF) 
S2(t)  which  was  (similarly  to  componCTits  chosen  previously  from  the  SDF  Sj^(t)) 
new  at  time  zero. 

The  process  continues  until  each  position  is  assigned  a different 
servicing  component  - corresponding  to  sub- sets  - and  containing  a single 


individual. 
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Vfe  are  interested  in  (i)  the  times  at  which  original  positions  receive 
a component  not  shared  with  any  other  position  and  (ii)  the  failure  times 
of  these  components. 

The  model  we  have  described  above  might  be  used  to  represent  the  failure 
times  of  elements  of  an  organism  in  vifhich  replacement  of  damaged  elements 
becomes  progressively  more  specialized.  The  first  failure  of  an  individual 
in  a sub^(single  individual) -set  might  be  regarded  as  corresponding  to  the 
failure  of  the  vAiole  system. 

System  (Procedure)  2:  Tliis  system  is  one  possible  inverse  of  System  (Procedure) 
1.  ’.'/e  suppose  that  initially  (at  time  zero)  we  have  k elements  (chosen  at 

I**  1 

random  frc«n  a population  with  SDF  S^(t))  grouped  in  sub  - sets  (as  in  the 

r*l 

penultimate  stage  of  System  (Procedure)  1) . V^hen  any  element  in  a sub  - set 
fails,  all  elements  in  that  sub^-^set  are  replaced,  in  effect,  by  one  servi- 
cing component  chosen  at  random  from  a population  with  SDF 
call  such  a component  a second  stage  component.  The  sub  - sets  are  grouped 
in  sub^-'^sets.  A sub  - set  may  contain  some  sub  sets  sendced  by  a coranon 
(second  stage)  crarqponent,  and  some  in  which  the  original  elements  still  survive 
(not  having  failed  as  yet).  Should  any  of  the  second  stage  components  fail,  all 

elements  in  tlie  sub^-^set  to  which  it  belongs  are  replaced,  in  effect,  by  a 

third  stage  conpment  servicing  every  position  in  the  sub  - set.  The  sub  - sets 

are  grouped  in  sub^-  sets,  and  so  on.  Generally,  a sub  - set  can  contain 

sub^’-^^sets  each  having  a common  servicing  component,  for  j=l,2,...,s  (a 
sub-set  being  * as  indicated  above  - a single  elonent) . VJlien  any  corponent 
servicing,  a sub’"'-'*’^set  fails,  all  positions  in  the  sub^-^set  to  whidi  it 
belongs  are  serviced  by  a common  ((s+1)- stage)  component,  chosen  from  a 
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population  with  SDF  S^_g(t). 

As  in  System  (Procedure)  1,  all  corponents  are  new  at  time  zero. 

Eventually  we  should  reach  a situation  where  v/e  have  some  (r-l)-th  stage 
elements  present.  These  are  not  members  of  any  subset,  and  when  one  of  them 
fails,  all  elements  (of  all  stages)  have  to  be  replaced. 

This  model  corresponds  to  progressively  less  specialized  replacement. 

The  time  of  failure  of  the  first  (r-l)-th  stage  element  to  fail  might  be 
regarded  at  the  time  of  failure  of  the  organism  as  a whole, 

3.  Notation 

V7e  now  develop  a notation  idiich  will  enable  us  to  express  the  operation 
of  Procedures  (Systems)  1 and  2 more  precisely.  The  formalization  devised 
below  can  be  applied  with  equal  facility  to  some  other  (obvious)  variants  of 
the  two  extreme  models  described  in  the  previous  Section.  An  example  v;ill  be 
given  at  the  end  of  this  Section. 

It  is  first  necessary  to  establish  identification  of  various  sets  of 
elements.  Vfe  siropose  that  the  k positions  are  initially  split  up  into  s^^ 
"first-level”  (sub)sets;  tliat  the  i-th  first-level  set  is  split  up  into 
S2(i)  second- level  (sub  )sets;  that  the  j-th  second  level  set  in  the  i-th 
first  level  set  is  split  up  into  S2(i,j)  third  level  (sub ^) sets,  and  so  on. 
Generally,  we  denote  the  ij^-th  h-th  level  set  in  the  ifj.i’th  (h-1)  level 
set  in  ...  in  the  i^-th  first  level  set  by  (ij^,i2, . . . ,ivj)  for 
ij^  * I2  “ l,...,S2(rj^)j  i^  *"  and  " 1 , . . . , 

®h  ^^1*^2’ ■ ■ ■ *^h-l^ ’ where  h runs  from  1 to  r. 

Figure  1 sets  out  the  corresponding  subdivision  diagrararkatically.  Note 
that  in  System  (Procedure)  1 the  "level"  is  the  same  as  the  "stage"  but  in 
System  (Procedure)  2 the  relation  between  stapes  and  levels  is  given  by: 
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stage  = r+l  - (level) . 


Finally,  we  denote  the  number  of  elanents  in  Sj^  (ij^,i2,. . . ,ij^)  by 
kh(ii,i2»- • Clearly  the  sum  of  kjjCO  over  all  combinations  of  values 
of  the  arguments  is  equal  to  k,  (the  number  of  the  positions)  that  is 

1 I * * ■ *^h^  = k . ( 

h "h 

Also  it  folloi\rs  from  the  definition  that 


’■■'hole 


sub-sets 


sub-sets 


S^d) 

I 


5^(2) 


S^(si) 


(1 ,1)  • • • ■^  (1  ,S2  (1))  ...  ^ (Sj  ,1)  • • »S2  (s2^  ,S2  (Sj) ) 


. .”;n 

S3(A2,1) .. .S3(A3) 


sub-sets  S^(l,l,l) .. .52(1,1,53(1,1))  . 


sub-sets  5 (1,1,. . .1) . . . 5 (1,1,. . . 
(elements)  . . .s^(l,. . .1)) 


[~TT":  I 

5j.(A^.l,l)...S^(Ap 


FIGURE  1.  Sdiematic  representation  of  levels  and  stages  ir.  hierarchal  replacement 
models. 

Note:  A.  is  defined  by  t)ie  recurrence  relation 


with  • 
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i: 


Using  this  notation,  System  (Procedure)  1 can  be  defined  succinctly  as 
follo\/s: 

(i)  Initially  (at  time  zero)  we  have  a set  Sq  of  Ic  elements , one  in 
each  of  k positions,  all  serviced  by  a single  corponent  randomly  chosen  frm 
a population  with  survival  distribution  function  (SDF)  SQ(t); 

(ii)  Ii/hen  the  servicing  component  for  a sub^set  Sj('>)  (j=l,...,r)  fails, 
it  is  replaced  by  Sj^j^('’)  different  servicing  components,  one  for  each  of  the 
sub^-^sets  in  Sj(»)-  Tliese  servicing  components  are  chosen  at  random  from 
the  survivors  in  a population  with  SDF  Sj^^(t)  which  was  new  (i.e.  originated) 
at  time  zero. 

System  (Procedure)  2 can  nov/  be  defined  as  follows: 

(i)  Initially  (at  time  zero)  we  have  1:  individual  elements  ("siib-sets") 
each  c’losen  at  randcm  from  a pooulation  with  SDF  Sj.(t) . 

(ii)  ’/.hen  the  servicing  conponent  for  a sub-set  ~j(0  Cj=r,  r-l,...,l) 
fails,  the  sub^-^aet  ■Sj_j^(°)  to  v/hich  it  belongs  is  assigned  to  a cannon 
servicing  component  chosen  at  random  from  the  survivors  of  a population  with 
SDF  Sj_j^(t)  Avhich  was  new  at  time  zero.  This  component  now  replaces  all 
corponents  previously  servicing  elements  belonging  to  5j_j^(°)* 

"h  may  have  a irode'l  in  vhich  the  two  systems  are  mixed.  A 
siiple  case  would  be  vhen,  initially,  replacement  is  according  to  System  1, 

but  after  reaching  the  stage  of  stib^sets  with  separate  servicing  ccmnonents 

(h  being  specified)  a failure  of  any  of  these  components  results  in  assignment 

’i-l 

of  a newly  selected  servicing  component  to  the  \diole  of  the  svh*  - set  to 

U 

which  the  corresponding  sub- set  belongs  (as  in  System  2).  Subsequently, 

System  1 is  followed  throughout  (so  that  the  next  time  a sub-set  servicing 


l! 


1 1 


t I 


i 

t 
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component  fails  each  of  the  sub^-^sets  contained  in  the  sub-set  is  assigned 
a different  freshly  chosen  servicing  compcnent). 

We  could  denote  the  sets  involved  in  the  first  period  (v/hen  System  1 
is  operating)  by  5^(0.  and  those  thereafter,  by  SjCO>  starting,  of  course, 
mth 

Further  generalizations  to  a number  of  changes  from  System  1 to  System  2 
and  conversely,  could  be  acconrwdated  by  including  greater  detail  in  the 
superfix  of  S("),  or  by  using  prefixes.  For  exairple  if  there  were  changes 
from  System  1 to  2 on  failure  of  h2^,h2,. . .,h^^  level  corponents,  and  converse 
changes  at  22,g2,...,gq  with 

Cl  < hi  < g2  < V-*'%  \ 

(no  cliange  in  system  at  second  failure  at  an  h^^  stage)  the  set  Sj , for 

< j hj^  readied  after  the  change  (from  System  2 to  1)  at  the  Ci'th 

level  v/ould  be  denoted  by  S.(")- 

®i  ^ 


4.  Distributions  of  Failure  Times;  System  (Procedure)  I . 

"fe  first  consider  the  failure  time  distribution  for  the  servicing 
coirponent  of  the  sub-set  Sj  (1,1,...,  1).  This  failure  tine  wall  be  denoted 
by  Tj.  (If  necessary,  the  more  explicit  notation  Tj  (1, !,...,!)  could  be 
used.)  The  sub^set  Sj  (1,1,...,!)  contains 


I ”*■’  I S (!,...,!, a.^j,...ap 


(4) 


components . 

Let  (Yj)  be  a sequence  of  independent  random  variables  (j=0,l,...)  wdth 
SDF’s  S^(t). 

’/'hen  Sq  fails  at  time  yp  (according  to  EDF  SQ(yQ)),  the  replacement 
servicing  compment  (for  the  subset  Si(l)  containing  S-  (1,1,...,!))  has  a 


s 


conditional  failure  time  distribution  as  that  of  random  variable  Y,  truncated 
from  below  at  Xq.  The  conditional  5DF  is 

{SiCXq)  'Sj(t)}/S^(yo)  (t^yg)  C5) 


and  the  corresponding  conditional  density  is 

flCt)/S^(yQ)  (tsyg) 


(6) 


where  fj(t)  = -dSj(t)/dt  . 

Hence  (with  yg,yj^, . . . ,yj_j^,yj  denoting  failure  times , at  successive  stages, 
of  the  servicing  conponents,  for  the  Whole  set,  subset, ... ,sub^-^set  and 
sub^set  (which  is  just  Sj(l,l,. .. ,1))  respectively  containing  the  subset 
Sj  (1,1,...,!))  we  have 


T ^j  ^2  ^1  f,(yT)  f.(yj 

Pr[Ti^T]  = / / •••  / / fv  1 s fy^ 

1 0 0 0 0 ^^  ^i^V  j'-^j-i^  ^ 

= dy  dy  . . .dy  (7) 

0 0 0 0 n?,jS.(y,_p  " 1 } 

T ^2  ^1  j-1  ^4(74) 

= / / ...  / / {n  g ^}f.(yjdyg<h/^...dy.  . (7)' 

0 0 0 0 i=0  ^i+l'-^i-'  ^ ^ ^ J 


For  j = 1 we  have 

T >^1 

Pr[Tj^ST]  = / fj(y^)  / {fQ(yg)/S3^(yQ)}dyQdyj^ 

= / ^fo‘^yo^/^i^>'o^^  / fi(yi)dyidyo 

0 yg 

“ ^ tfo^^o^/^l^yo^^^^l^yo^  ‘ S2(T)}dyQ 


= 1 - Sg(T)  - Sj^(t)  nfo^yo^/^i^yo^^^yo 
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Pr[T,>T]  = SqCt)  - 5,(1)  / 

X U X Q 


T dS^Cyg) 


in  agreement  with  Krakowski's  formula  (1.1). 

It  is  interesting  to  note  that  the  distribution  given  by  (8)  differs  from 
the  distribution  of  age  at  failure  of  a randcmly  chosen  individual  from  Sj^(0» 
given  that  it  exceeds  age  at  failure  of  a randomly  chosen  individual  from 
Sr,(“).  The  cdf  of  the  latter  is 


/ f,(t)FQ(t)dt 
0 

^^55 , (8a 

/ £i(t)Fo(t)dt 

where  F^(t)  = 1 - Sj^(t),  i=0,l. 

In  (8)  we  suppose  we  obtain  an  observation  corresponding  to  eoery  age 
at  death  of  tlie  indiAridual  from  Sq(->);  in  (8a)  an  observation  is  obtained 
only  when  the  individual  from  Sj^(»)  survives  that  from  Sq(0,  so  that  the 
larger  values  of  the  argument  get  relatively  low  lighting.  Note  also  that 
failure  times  for  any  individual  servicing  conponent  (including  that  of  the 
elements  in  the  final  stage)  are  not  independent;  therefore  the  distribution 
of  tte  first  or  the  last  failure  time  among  these  components  or  elements  is 
not  easily  obtained. 


Hxaroles.  a)  If  we  suppose  that  each  servicing  component  (at  whatever 
level)  has  a Weibull  distribution,  with  SDF 

Sj(t)  = exp{-(t/e)‘^}  (UO;  9 >0)  (9) 


so  that 
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-li  >^2^1  j 

■U  j , , J , * , 


then  (7)  becomes 

Pr[Tj^T]  = (cO  / / .../  / ( n (y  /e)}^‘-‘^exp{-(yye)‘^}dy„dy,..  .dy. 

0 0 0 0 i=0  ^ ^ y i 1 

= 0C- V' . . ./'V' ('i'ti)  . . . dt 

0 0 0 0 i=0  ^ ^ ^ 3 

= -1 

j ! 0 •'  ^ 

In  the  case  given  by  formula  (11)  T^/e  has  a generalized  gamma  distribu- 
tion (see,  e.g.,  Johnson  and  Kotz  (1970)).  It  is  distributed  as  the  c-th 
root  of  a gamma  variable  vdth  parameter  (j+1) , 

b)  Taking  c = 1 gives  the  expression  for  the  special  case  of  coninon 
exponential  survival  distribution  functions. 

c)  If  we  suppose  tliat 


-t/e 


S^(t)  « e ^ (tsO;  e^>0) 


(12) 


with  not  all  6 s equal,  then  we  find 
J .1  T >^j  >^2  ^1 


Pr[T.sT]  = ( n / •../  / exp{- J (e:^-6:J;py.-e:V,-}dyodyi...dyi> 

■’  1=0  ^ 0 0 0 0 i=0  ^ ^ ^ ^ J J ^ ^ J 


Integrals  of  the  form 


T yj  '^2  ^1  1 

I (an. a, a.)  = / / .../  / exp(-  \ a.y.)dyndy, . . .dy. 

^ 00  00  i^o^’-  ^ J 


0 0 0 0 

satisfy  tlie  recurrence  relation: 


(13) 


(14) 


1 » • • • 9 aj ) * aQ  {I  (a^^ >^2  > • * • ” 1 ^ a^'^a^ , a^ , • • • >a ^ ) } . (15) 


Using 


T 1 -a,.T 

I(aQ)  = / exp(-aQyQ)dyQ  = a‘^(l-e  ) 


(16) 
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[ 


and  repeatedly  applying  (15)  vre  can.  successively  obtain  the  explicit  values 
of  I C^Q  I ^2  > • • • > j ) • 

Introducing  the  displaceinent  operator  E,  such  t!iat  E^f(x)  = f(x+h)  we 


have 


J .1  J a. 

ICao,a, ,...,a.)  = ( n a.)  ^ n (1-E  ^)exD(-O.T) 

^ ^ J i=0  ^ i=0 

i 1 j a-  a.+a., 


i'O 


i=0  i<i' 


+ (-l)^‘^^E^^^^]exp(-0.T) 
-(a-+a. ,)t 


J -1  1 ■a.j'r  -(a.+a- 

= C n aj  ^{1-  I e ^ + n e " ' 
i=0  ^ i=0  i<i'  . 


i+1 

+ (-1)^  ^e  ^ ^ . (17) 


In  System  (Procedure)  2,  the  situation  is  somewhat  more  coup  lex.  H'W- 
ever,  the  distributions  o£  final  failure  time  can  be  explicitly  derived  without 
much  difficulty.  For  coirputational  sinplicity,  we  will  restrict  our  calcula- 
tions to  the  symmetric,  balanced  case,  in  which  the  number  of  sub^sets  in  each 
sub^-^set  (i.e.  serviced  by  the  same  (j-l)-th  level  conponent)  is  the  same  - 
Sj  - say.  In  this  case  k = and  the  number  of  sub-sets  is 

We  denote  the  survival  time  distribution  function  for  an  h-th  level 
service  conponent  (failing  in  service)  by  (t) . This  is  to  be  distinguished 
from  the  SDF  Sj^(t)  in  the  population  from  which  the  h-th  level  service  com- 
ponents are  chosen  when  a relevant  (h+l)-t.h  level  corponent  fails.  From  the 
syranetry  of  the  hierarchal  structure,  the  SDF  for  each  of  the  s^^  h-th  level 
corponents  to  be  serviced  by  any  one  (h-1)  level  component  is  the  same,  and 
so  the  SDF  of  the  first  failure  time  among  these  Sj^  h-th  level  components  is 
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Conditional  on  this  time  beint?  equal  to  yjl,,  say,  the  SDF  of  failure  time  of 
the  (h-l)-tl\  level  component,  activated  at  time  y,'  is  the  truncated  distri- 


bution 


SivlW 


(t^yA)  . 


Averaging  over  the  distribution  (19)  of  first  failure  time  of  the  relevant 
h-th  level  conponents,  we  obtain 

t {s.--^(y)}"^  ^fp.^(y) 

2(h-i)  w = WiW/  - 


vhere  fr-^y)  = - dS^s(y)/dy  . 


Vq  also  have  for  the  ultimate  level  r: 


S(r)(t)  = Sr(t)  . (22 

Applying  (22)  repeatedly  with  h = r,  r-l,...,l  we  eventually  obtain  an 
explicit  formula  for  S^g^(t),  the  SDF  of  'final  failure  time." 

y\l though  this  calculation  can  be  effected  numerically  in  a straight- 
forward v/ay,  it  often  leads  to  rather  conplicated  algebraic  ejq^ressions . 

In  the  follo\'ring.  derive  an  expression  for  ^‘’(•^.-2)^^^  where  each  SDF  is 
a 1‘feibull  with  the  same  c,  but  different  9 parameters. 


Tnen,  from  (21)  and  (22) 


(ex,  -W 


t [exp{-(y/ej}^]  ce'^(y/6_)^'^exp{-(y/0j^} 


exn{-(y/e^_j^)  } 


=\exp{-(t/e^_^)^}/  c9‘  (y/0j.)  ■ exp[-{Sj.e'  - 6’.;i}y  ]dy 
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VfVrV 

If  = (e^9^.p 


j— ^-[exp{-(t/e^_^)  } - exp{'Sj.(t/0j.)  }] 


S(r-i)Cf)  “ . 


Again  using  (21)  and  (22) , we  find  (assuring 


^(r-2)^^^ 


s s , 

r 1 r*l 


' s^.lOxp(-Ct/(1^.2)  1 • ( 1-1  cl 

Sr-CVr-d 

f s , -1 

X / [exp{-(y/0^_j)^}  - exp{-Sj.(y/0p^}] 

cy^’^[9*?fxp{-(y/6r,i)'^}  - s^0‘%xp{-s^(y/op'^}] 

exD{-(y/0^_2)'^J 

' ^r-lirT/eT'-cl  '■■W-(t/6^.2)''} 

Sr  tOr0r-i3 

t ^ i^r-l"^  . s T-1 

.c-1  V /•  1^1  r ^"i  ■> 


(23.1) 


(23.2) 


i[  ^ )[9;fiexp{-((s^.^-i)0;^--js^e;''-0;f2)y''} 


- s^o;"exi'.{-((s^.j-j-i)o;^;^  - (j+Ds^©;""  - 0;?2)y‘'>]dy 

f iVAy^  \ nJrVr^r  -c 

^r’(Vr-l^  3-0  I ^j‘®r-2 


- s 0 
r r 


Bi*re;-2  ’ 

lAere  = Cs^_j-j)0^^j  - 35^.0'*^  provided  f 0*^2  ^y  i 


(24) 


3 


is  replaced  by 

tSxp{-(t/0^_2)^}. 


The  ^'non-symmetria"  case  reqvdres  somewhat  more  involved  comnutational 
and  teclmical  analysis. 

Hie  relevant  formula  for  the  survival  distribution  of  the  cc»nponent 
servicing  Vl^^l’*“’Vl^  “ 


(t)  = S.  ,(t)/  1- 
n I Q 


^lufh) 

aTh-i) 


(25) 


where  uCh)ei(h-l)  means  that  uGO  = (i]^,i2»- • • :.ij^.]^»UjP  • 
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